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(^ Abstract 

CnJ Common intervals of K permutations over the same set of n elements were firstly investigated by 

Vh T. Uno and M.Yagiura (Algorithmica, 26:290:309, 2000), who proposed an efficient algorithm to find 

Mh common intervals when K — 2. Several particular classes of intervals have been defined since then, 

■^ e.g. conserved intervals and nested common intervals, with applications mainly in genome comparison. 

(yQ Each such class, including common intervals, led to the development of a specific algorithmic approach 

T— I for K — 2, and - except for nested common intervals - for its extension to an arbitrary K. 

In this paper, we propose a common and efficient algorithmic framework for finding different types 
of common intervals in a set V of K permutations, with arbitrary K. Our generic algorithm is based on 
a global representation of the information stored in V, called the M/nMax-profile of V, and an efficient 
data structure, called an LR-stack, that we introduce here. We show that common intervals (and their 
subclass of irreducible common intervals), nested common intervals (and their subclass of maximal 
nested common intervals) as well as conserved intervals (and their subclass of irreducible conserved 
intervals) may be obtained by appropriately setting the parameters of our algorithm in each case. All the 
resulting algorithms run in 0{Kn + A^)-time and need 0{n) additional space, where N is the number 
^•^ of solutions. The algorithms for nested common intervals and maximal nested common intervals are 

T^lj- new for K > 2, in the sense that no other algorithm has been given so far to solve the problem with the 

T-H same complexity, or better. The other algorithms are as efficient as the best known algorithms. 
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1 Introduction 

Common, conserved and nested common intervals of two or more permutations have been defined and 
studied in the context of genome comparison. Under the assumption that a set of genes occurring in 
neighboring locations within several genomes represent functionally related genes ifTTl [1711211 . common 
intervals and their subclasses are used to represent such conserved regions, thus helping for instance to 
detect clusters of functionally related genes |:18 , 22J , to compute similarity measures between genomes 
161 or to predict protein functions |[T5ll24l . Further motivations and details may be found in the papers 
introducing these intervals, that we cite below. 

In these applications, genomes may be represented either as permutations, when they do not contain 
duplicated genes, or as sequences. In sequences, duplicated genes usually play similar roles and yield the 
interval search more complex |fT0ll20l , but sometimes they are appropriately matched and renumbered 
so as to obtain permutations [91 12] . 

We focus here on the case of permutations. Efficient algorithms exist for finding common and con- 
served intervals in K permutations (K > 2), as well as for finding irreducible common and irreducible 
conserved intervals |23 , 12 , 5,61. Nested common intervals and maximal nested common intervals have 
been studied more recently lfT4l . and efficient algorithms exist only for the case of two permutations ID. 
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Surprisingly enough, whereas all these classes are subclasses of common intervals, each of them 
has generated a different analysis, and a different approach to obtain search algorithms. Among these 
approaches, interval generators (5] have been shown to extend from common intervals to conserved 
intervals |19|, but this extension is not easily generalizable to other subclasses of common intervals. 

The approach we present in this paper exploits the natural idea that an efficient algorithm for com- 
mon intervals should possibly be turned into an efficient algorithm for a subclass of common intervals 
by conveniently setting some parameters so as to filter the members of the subclass among all common 
intervals. It also chooses a different viewpoint with respect to the information to be considered. In- 
stead of searching intervals directly in the permutations, it first extracts the helpful information from 
the permutations, focusing on each pair (i, i + 1) of successive values in {1, 2, . . . , n} and defining 
the so-called M/nMax-profile of the permutations. Then, it progressively computes the set of interval 
candidates, but outputs them only after a filtering procedure selects the suitable ones. 

The organization of the paper is as follows. In Section |2] we present the main definitions and the 
problem statement. Then, in SectionlS] we introduce the abstract data structure on which strongly relies 
our main algorithm called LR-Search, also described in this section. The complexity issues are discussed 
in Section!?] In Section l5] we give the specific settings of our algorithm for common, nested common 
and conserved intervals, and prove the correction in each case. In Section l6] we show how to further 
modify the algorithm so as to deal with even smaller subclasses. SectionJTjis the conclusion. 

2 Generalities 

For each positive integer u, let [u] :— {1, 2, ... , u}. Let V := {Pi, P2, ■ ■ • , ^a'} be a set of permuta- 
tions over [n], with n > and integer. The i-th element of P^ is denoted p*^, for all i E [n]. The interval 
[i,j] of Pk, defined only for 1 < i < j < n, is the set {Pi,Pi^i, ■ ■ ■ ,p'j}, and it is denoted {i..j) when 
Pfc is the identity permutation /c/„ := (12 ... n). Then {i..j) = {i,i + 1, . . . ,j}. For an interval [i,j] 
of Pfe, we also say that Pfc is delimited by p^ and p'^, or equivalently that p^ and p^' are the delimiters 
of [i,j] on Pfc. Furthermore, let P/^^ : [n] — > [n] be the function, easily computable in linear time, that 
associates with every element of Pk its position in P^,. We now define common intervals: 

Definition 1. f23l A common interval of T' is a set of integers that is an interval of each Pk, k £ [K]. 

Nested common intervals are then defined as follows: 

Definition 2. lfT4ll A nested common interval (or nested interval for short) of T' is a common interval / 
of V that either satisfies |/| = 2, or contains a nested interval of cardinality |/| — 1. 

Example 1. Let Pi = W7 and P2 = (7213645). Then the common intervals of "P = {Pi, P2} are 
(1. . 2), (1. . 3), (1. . 6), (1. . 7), (3. .6), (4.. 5) and (4..6), whereas its nested intervals are (1.. 2), (1.. 3), (4.. 5), 
(4..6)and(3..6). 

With the aim of introducing conserved intervals, define now a signed permutation as a permutation 
P associated with a boolean vector signp that provides a + or — sign for every element of P. Then 
signp [i] is the sign of the integer i in P. Note that, even in a signed permutation, the elements are 
positive integers. An element of P is called positive or negative if its associated sign is respectively 
+ or — . A permutation is then a signed permutation containing only positive elements. Conserved 
intervals need to assume that p\ — 1 (positive) and p^^ n (positive), for each k G [K]. 

Definition 3. [|61| Let V := {Pi, P2, . . . , Pk} be a set of signed permutations over [n], with p'^ — 1 
(positive) and p^ — n (positive), for each k E [K] . A conserved interval of T-" is a common interval of 
■p (ignoring the signs) delimited on Pk by the values a^ (left) and bk (right), for each k e [K], such that 
exactly one of the following conditions holds for each k: 

(1) ttk ~ fli with the same sign and bk = &i with the same sign. 

(2) ttk — bi with different signs and bk = ai with different signs. 



Example 2. Let Pi = Id-; and P2 = (13 2 6457) with signp^ = [+,-,-,-,-,+,+], mean- 
ing that 1,6 and 7 are positive, whereas 2,3,4 and 5 are negative in P2. Then the conserved in- 
tervals of T' = {Pi,P2} are (1..7) and (2. .3), whereas its common intervals (ignoring signs) are 
(1..3), (1..6), (1..7), (2..3), (2..6), (2..7), (4..5), (4..6), (4..7). 

In the following problem, C refers to a class of common intervals, e.g. common, nested or conserved. 

C-lNTERVAL Searching Problem (abbreviated C-ISP) 

Input: A set T' = {Pi , P2, ■ ■ ■ , Pk} of signed permutations over [n], satisfying the conditions required 

by the definition of C. 

Requires: Give an efficient algorithm to output without redundancy all C-intervals of P. 

In all cases, we may assume without loss of generality that Pi = Id^, by appropriately renumbering 
Pi . The other permutations must be renumbered accordingly. 

In this paper, we propose a common efficient algorithm to solve C-ISP when C stands for common, 
nested, conserved intervals as well as for their respective subclasses of irreducible common, maximal 
nested and irreducible conserved intervals (defined in Section[6]l. For common and conserved intervals 
(irreducible or not), efficient algorithms have been proposed so far, developping different and sometimes 
very complex approaches. For nested and maximal nested intervals, efficient algorithms exist for the 
case K = 2. We solve here the case of an arbitrary K. Our approach is common for all classes, up to 
the filtering of the intervals in C among all common intervals. 

3 Main Algorithm 

Our LR-Seaich algorithm is based on two main ideas. Fkst, it gathers information from V that it stores 
as anonymous constraints on each pair (t, t + 1) of successive values, since during the algorithm it is 
useless to know the source of each constraint. Second, it fills in a data structure that allows us to find 
all common intervals with provided constraints if we need, but an additional Filter procedure is called 
to choose and to output only the intervals in the precise class C for which the algorithm is designed. 

The algorithm uses a specific data structure that we call an LR-stack. The candidates for the left 
(respectively right) endpoint of a common interval are stored in the L (respectively R) part of the LR- 
stack. At each step of the algorithm, the LR-stack is updated, the solutions just found are output, and 
the L/?-stack is passed down to the next step. 

3.1 The L7?-stack 

The L/?-stack is an abstract data structure, whose efficient implementation is discussed in Section]?] 
Definition 4. An LR-stack for an ordered set S is a 5-tuple {L, R, SL, SR, R^) such that: 

• L,R are stacks, each of them containing distinct elements from E in either increasing or decreas- 
ing order (from top to bottom). The first element of a stack is its top, the last one is its bottom. 

• SL, SR C 'S respectively represent the set of elements on L and R. 

• R^ : SL -^ SR is an injection that associates with each a from SL a pointer to an element on R 
such that R^ (a) is before R^ {a') on R iff a is before a' on L. 

According to the increasing (notation +) or decreasing (notation -) order of the elements on L and 
R from top to bottom, an L/?-stack may be of one of the four types L^R^,L^R^ , L^R^ , L^R^. 

Remark 1. We assume that each of the stacks L, R admits the classical operations pop, push, and that 
their elements may be read without removing them. In particular, the function topQ returns the first 
element of the stack, without removing it, and the function next{u) returns the element immediately 
following u on the stack containing u, if such an element exists. 
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Figure 1 : M/nMax-profile of T' = {Id^ ,(721364 5)} with basic bounding functions b{t) = rrit and 
B{t) ~ Mt, and execution of the L7?-Search algorithm. The stack is initially empty. For each pair 
(i, ^ + 1) with i = 71 — 1, n — 2, . . . , 1, the corresponding L7?-stack (below the pair) is obtained from 
the preceding L7?-stack (on its right) by performing the sequence of operations written below the arrow 
linking the two Li?-stacks. 

We further denote, for each a E SL and with a' = next{a), assuming that next{a) exists: 

Setii{a) = {5 e 5/? I 5 is located oni? between/?^ (a) included and 7?^ (a') excluded} 

When next (a) does not exist. Self; (a) contains all elements between R^ {a) included and the bottom of 
R included. Then R^ (a) is the first (i.e. closest to the top) element of Setji{a) on R. 

The following operations on the LR-stack are particularly useful. Note that they do not affect the 
properties of an L/?-stack. Sets SetjiQ are assumed to be updated without further specification whenever 
the pointers R^ () change. Say that a' is L-blocking for a, with a' 7^ a, if a cannot be pushed on L when 
a' is already on L (because of the increasing/decreasing order of elements on L), and similarly for R. 

• Pop^{a), for some a e S: pop successively from L all elements that are L-blocking for a, push 
a on L iff at least one L-blocking element has been found and a is not already on L, and define 
R^ {top{L)) as top{R). At the end, either a is not on L and no L-blocking element exists for a, or 
a is on the top of L and R^ (a) is a pointer to the top of R. 

• Popj^{b), for some 6 e E: pop successively from R all elements that are _R-blocking for b, update 
all pointers R^ {) and successively pop from L all the elements a with R^ {a) = nil. At the end, 
either b is not on R and no i?-blocking element exists for b, or b is on the top of R. 

• Push]^fi{a, b), for some a, 6 e S (performed when no L-blocking element exists for a and no 
i?-blocking element exists for b): push a on L iff a is not already on the top of L, push 6 on i? iff 
b is not already on the top of R, and let R^ {top(L)) be defined as top{R). 

• FindL{b), for some 6 G S: return the element a of SL such that b E Set^{a). 

Example 3. In Figure [T| consider the L^i?+-stack below the pair (4,5). Here, Set^iA) — {5,6} 
and Seti^{l) ~ {7}. The instruction Popj^{3) discards 4 from L, and pushes 3 instead, also defining 
R^ (3) = 5. Next, the instruction fop^(6) discards 5 from R. The resulting LR-stuck is represented 
below the pair (3, 4). 



3.2 The Li?-Search algorithm 

Let 7^ be a set of K signed permutations over [n]. Recall that Pi = W„. In the following definitions, 
2<A;<iCandl<i<n-l: 



Algorithm 1 The LR-Search algorithm 



Input: Set V of signed permutations over [n], bounding functions b and B, Fiher procedure 
Output: All common intervals {t..x) ofV with x E Setu{t), filtered by Filter 

1: Compute nif, nit and bt with 2 < k < K and t E [n — 1] 

2: Compute Mf,Mt and Bt with 2 < k < K mdt e [n - I] 

3: Initialize an L^i?+-stack with empty stacks L, R 

4: for t from n — 1 to 1 do 

5 : Popi^ [bt ) II discard from L all candidates larger than bt and push bt instead 

6: Popfj {Bt ) // discard from R all candidates smaller than Bt 

1: ifBt=i+lthen 

8: Push]^ji{bt,t + 1) lit + 1 is a new right candidate, suitable for each a on L 

9: end if 

10: Call Filter to choose a subset of intervals {t..x) with x G Setji{t) 

1 1 : end for 



mf^ :— Ttmi{p^ \ pf is located between t (included) and t + 1 (included) on Pk} 

rrit := min{m^ \2 < k < K} 
Mt := max{p^ | p'^ is located between t (included) and t + 1 (included) on Pk} 
Mt -.^ max{M^ \2 < k < K} 
Notice that, for each k E [K], 

nit < m'l < t < t + 1 < M'I < Mf 

Now, let the bounding functions b,B : [n] -^ [a] be two functions such that b{t) < nit and B{t) > 
Mf Denote bt := b{t) and Bt := B{t) for all tE[n- 1]. 

Definition 5. Let P be a set of permutations on [n] . Then the MinMax-profile of V with respect to b 

and B is the set of pairs [bt^Bt], t E[n— 1]. 

The L/?-Search algorithm (See AlgorithmfTli works as follows. Each pair (t, i + 1), 1 < < < n — 1, 
is associated with the couple [fef,Bt], which means intuitively that every common interval {a..c) of T' 
that contains t and < + 1 should satisfy a < bt < Bt < c. The L~i?+-stack, initially empty, stores on 
L (respectively on R) the candidates for the left endpoint a (respectively right endpoint c) of a common 
interval {a..c), and links together the left and right candidates to form suitable pairs. Intuitively, a 
suitable pair is a pair of endpoints for which no disagreement has been found yet. In the LR-stack, a 
suitable pair is a pair (a, c) E SL x SR such that a < FindL{c). Equivalently, right candidate c could 
form common intervals with FindL (c) as well as with all left candidates following FindL (c) on L. 

Remark 2. It is important to notice that, by definition, in an L^i?+-stack, a' is L-blocking for a iff 
a' > a, and b' is i?-blocking for b iff b' < b. 

In its most general form, the algorithm outputs all common intervals without any redundancy. How- 
ever, in order to insure its efficiency even for strict subclasses C of common intervals, it has two types 
of parameters allowing us to filter intervals during the search (bounding functions achieve that) and 
during the output step (Filter procedure achieves that). Then, LR-Search considers (step 4) every pair 
(i, t + 1) using decreasing values of t, so as to output the intervals with same left endpoint altogether in 
the same step (increasing values of t would produce the intervals with same right endpoint in the same 
step). Operations Pop ]^{bt) , Pop j^{Bt) (steps 5, 6) successively allow us to discard left candidates a that 
do not satisfy a < bt (respectively right candidates c that do not satisfy Bt < c), and to update the 
suitable pairs accordingly. Steps 7-9 identify the cases where t + 1 may be added to R, namely when 



Bt — t + 1 {PushLn tests whether t + 1 is already on R). Finally, in step 10 suitable intervals {t..x) are 
necessarily common intervals since further constraints on pairs (t', t' + 1), t' < t, will not affect {t..x). 
The Filter procedure is then called to filter among all common intervals found here those that belong to 
some specific subclass. 

Example 4. Figure [T] shows the application of the L7?-Search algorithm on the set V ~ {Pi, P2} of 
permutations, where Pi = Idj and P2 = (7213645). In this case, we defined bt = nit — rnf and 
Bt = Mt = MJ. The constraints [bt^Bt] are written above the arrow representing the pair (t, t + 1), for 
alH < n — 1, so that the upper part of the figure represents the MinMax-piofile of V. The first step 
corresponds to t = 6, and consists in performing operations Popj^{l), Popj^{7) and PushLji{l,7) on 
the initially empty stack. Thus the first two operations have no effect, whereas the third one pushes 1 
on L,7 on R and defines R^ {1) as a pointer to 7. The next step takes the current state of the LR-stuck 
and performs Popj^ {A),Popj^{6) and PushLi^{4:, 6) to obtain the LR-stack below the pair (5, 6). The first 
common intervals are output when t — A, namely (4.. 5), (4.. 6). 

Remark 3. In the rest of the paper, the notation forj concerns the execution of the for loop in LR- 
Search for some fixed t. Similar notations will be used for the loops in the Filter procedures given 
subsequently. When these notations are not confusing, we use them without any further specification. 

Let 5ef^(a) be the value of Set]i{a) at the end of step 9 in forj, for each a on L. Let Pairst be the 
set of pairs {t..x) with x e Set]i{t). 

Theorem 1. Assuming the Filter procedure does not change the state of the LR-stack, the set A defined 
as A : — Ui<t<nPoirst computed by LR-Search is the set of all common intervals {t..x) ofV (ignoring 
the signs) satisfying t = bt and x = B^-i- Moreover, for all these intervals we have 

bt ~ niin{6tu \t < w < x ~ 1} 

Bx~i — niax{i?u, I i < ui < a; — 1}. 

We first prove the following result. Notations min{I) and max{I) for an interval I of P respectively 
denote the minimum and maximum value in /. 

Claim 1. Let t, 1 < t < n. After the execution of step 9 in for;, we have x G Setp{a) ifft, x and a 
satisfy the three conditions below: 

(a) a = min{6.u, \t<w<x~l} 

(b) X — Bx-i = max{i?tu \t < w < x ~ 1} 

(c) for each k, 2 < k < K, the interval Ik of Pk delimited by the leftmost and rightmost values 
among i, i + 1, . . . , a; satisfies min{Ik) > a and max{Ik) = x. 

Proof of Claim [l] We use induction on t, with decreasing values. Let t = n — 1 and notice that the 
L/?-stack is empty up to step 7 in for„_i. 

Proof of "=>:" The only operation that can affect the empty stack is PusliLnibn-i, n), performed only 
when i3„_i = n. Then Set^~^{bn-i) = {n}. Thus a = bn-i,x — n and affirmations (a),(&) 
are verified. Now, the interval Ik delimited by n — 1 and n satisfies the required conditions since 
min{Ik) — niji-i > bn-i = a and max{Ik) = n = x. 

Proof of "<^:" Let a, x and Ik be defined as in {a)-{c). Because of i = n — 1, together with x < n 
(since x is an element of [n]) and with t < w < x — l,wt deduce that x ~ n and w = n — 1. Thus 
a — &„_i. By affirmation (6), i?„_i ~ n. In conclusion, we have fo„_i — a and i?„_i = n thus the 
conditions in step 7 are fulfilled. Then a is pushed on L, n is pushed on R and we are done. 

Assume now the claim is true for i + 1, where i + 1 < n — 1, and let us prove it for t. For that, we 
denote a' := inm{bw \t + 1 < w < x ~ I}. 

Proof of "^:" With x e Set]^{a), three cases are possible: 



(i) X e Setj^ (a). In this case, during fort a is not discarded by Popj^{bt), thus (1) a < bt. 
Moreover, by the inductive hypothesis (affirmation (a) for t + 1), x e 5ef^^(a) impHes that 
(2) a = inin{feu, \t + 1 < w < x — 1} = a' . Then by (1) and (2): 

iiiin{6tu \t < w < X — 1} = mm{bt, a} = a 

and affirmation (a) holds for t. Furthermore, x is not discarded by Popj^{Bt), thus (3) x > Bt. 
Affirmation (6) follows using the inductive hypothesis. Now, let /^ be the interval given by 
affirmation (c) for i + 1 and Pk- Let I'^ be the interval of Pk delimited by t and t + 1. Defining 
Ik — I'^U I'l, we have that Ik is an interval of Pk, since both /^ and I'l contain t + 1. Also, 
/j, is delimited as required, as /^ and I'^ do. Moreover, using the hypothesis bt < nit < nit '^^'^ 
Bt > Mt > Mt, properties (l)-(3) above and the inductive hypothesis on /^ we obtain: 

mm{Ik) — inin{min{I'f.),min{I'^)} > iiiin{a',W(} > inin{a,OTj} > min{a,^t} — a 
max{Ik) = max{mflx(/^), mflx(/^')} = max{x, Mt} — x. 

(ii) there exists a" > a such that x e Set*j^^{a"). Then a" is necessarily discarded by Popj^{bt), thus 
a" > bt and Set*j^^{a") goes entirely into Set]f{bt)- Then, we deduce x e Set]^{bt) and therefore 
X e Set\i{a)r\Set*ji(bt), which is impossible unless (4) a = bt- Now, by the inductive hypothesis, 
X e Set*j^^{a") implies by affirmation (a) that a" = min{6iu | i + 1 < w < a; — 1} = a'. Thus, 
recalling that a < a" by the hypothesis of case (ii), we have (5) a' ~ a" > a. With (4) and (5) 
we deduce: 

min{6^„ \t<w<x— 1}= imn{bt, a"} = imn{bt, a'} = minja, a'} = a 

and affirmation (a) holds for t. Furthermore, x is not discarded by Popj^{Bt), thus (6) x > Bt- 
Affirmation (6) follows using the inductive hypothesis. As before, let /^ be the interval given by 
affirmation (c) for i + 1 and Pk- Let I'^ be the interval of Pk delimited by t and t + I. Defining 
Ik — I'k^ I'k, we have again that Ik is an interval of Pk delimited as required. Moreover, using 
the hypothesis bt < nit < ni^ and Bt > Mt > M\, properties (4)-(6) above and the inductive 
hypothesis on /^: 

min{Ik) = uim{min{I'i^),min(I'^)} > niin{a",mj} > unn{a,mt} > uim{a,bt} — a 
max{Ik) — max{mflx(/^), mflx(/^')} — max{x, Mt} = x. 

(iii) there is no a" such that x G 5ef^^(a"). In this case, x is added to R during forj, in step 8. Then 
X ~ t + 1 = Bt, a = bt and affirmations (a) and (b) clearly hold. Recall that, by definition, 
t + 1 < Mt < Mt < Bt thus Mt = t + 1 = X. Now, the interval h dehmited by t and f + 1 on 
Pk satisfies: 

min(Ik) > nit > bt — a 
max{Ik) =Mt=t+l = x 
The "^" part of the claim is proved. 

Proof of "■^:" Let a, x and Ik,'2, <k < K, be defined as in affirmations {a)-{c) in the claim. Consider 
the two following cases: 



(i) t < X ~1. Notice that a' , x and the intervals I'f., 2 < k < K , delimited on Pk by the leftmost and 
rightmost values between t + 1, . . . ,x satisfy affirmations (a)-(c) for t + 1, thus by the inductive 
hypothesis x £ Set*j^^{a'). Now, we show that in both cases occurring during fort, we have 
(7) Set*j^^{a') C Set]^{a). Indeed, when a' is not discarded by Pop^{bt), we necessarily have 
a' < bt and thus we also have by affirmation (a) that a = min{&t, a'} — a'. Property (7) follows. 
When a' is discarded, we necessarily have bt < a' and thus, by affirmation (a), we deduce 
a = min{6f,a'} = bt- The execution of Pop j^{bt) implies, in this case, that Set*j^^{a') becomes 
a part of Set]^{a), and (7) is proved. Now, with (7) and given that, in step 6, x — B^-i > Bt 
implies that x is not discarded by Popj^[Bt), we obtain that x E Set'j^{a). 

(ii) t = X — 1. Then by hypothesis a = bx-i,x = B.j._i,min{If;) > b^-i and max{Ik) = x. In 
step 5 of for^;-!, the instruction Pop j^{hx-i) discards all a" > b^-i (if any). Thus, at the end of 
step 5, top{L) < bx-i- The instruction Popj^{x) insures that, at the end of step 6, top{R) > x. 
(Notice that x = t + 1 and all the elements pushed before on R by PusIilr are of the form t' + 1 
with t' > t). Then, the instruction PushLji{bx-i, x) pushes b^-i on L if necessary, pushes x on 
R (thus top{L) — bx-i and top{R) — x) and adds x to Set*'j^{bx-i). 

Claimlllis now proved. I 

Proof of Theorem^ By definition, {t..x) e Pairst iff a; G Set^j^{t). According to Claim [I] this holds 
iff affirmations (a)-(c) hold with a — t, that is, the following affirmations hold simultaneously; 

(a') t — min{6i„ \t < w < x — 1} 

(b') X = Bx-i = max{i?^ \t < w < x ~ 1} 

(c') for each k, 2 < k < K, the interval Ik of P delimited by the leftmost and rightmost values 
between t,t + 1, . . . ,x satisfies min{Ik) > t and max{Ik) = x. 

We show that Ik has the same elements as {t..x) and that bt = t (Bx-i = a: is directly given by 
affirmation {b')). By affirmation (c'), every element in {t..x) is an element of 7^. Conversely, assume 
by contradiction that u is any element of I^ distinct from t,t + 1, . . . , a;. Then u is not a delimiter 
of Ik- Consequently, let t', t < t' < x — 1, such that u is between t' and t' + 1 on Pk- Then 
bf < Jnt' < u < Alt' < Bf'. By affirmations (a') and (6'), t < bf and Bf < B^-i — x. Thus 
t < u < X, a contradiction. 

To show that bt — t, notice that bt < t, by the definition of bt, and t < bthy affirmation (a'). I 

4 Complexity issues 

We separately discuss the implementation of an LR-stack, and the running time of the LR-Seaich algo- 
rithm. 

4.1 The Li?-stack 

The efficient implementation of an LR-stack depends on the need (or not) to implement FindL- If Find l 
is not needed, then L and R may be implemented as lists. Consequently, Pop^ and Popj^ are easily 
implemented in linear time with respect to the number of elements removed respectively from L and R, 
whereas Push^ji takes constant time. Also, top{R),top{L) and nextQ need 0(1) time. 

When FindL is needed, then we are in the context of a Union-Find-Delete structure, where the 
operations are performed on the sets Setii{a), as follows: unions are performed by Popj^ and PusIilh, 
whereas deletions are performed by Popj^. These algorithms are already very efficient in the most 
general case [1|, but unfortunately not linear. Yet, particular linear cases may be found and show useful 
(see Algorithm [2|. 



4.2 The Li?-algorithm 

We prove the following result. 

Theorem 2. Assume that computing bt and Bt takes negligible time and space with respect to the com- 
putation of m^,mt,Mf and Mt insteps 1-2. Then the LR-Search algorithm has running time 0{Kn+ F) 
and uses 0{n + f) additional space, where F and f respectively denote the running time and additional 
space needed by the Filter procedure, over all values oft. 

Proof of Theorem l2l Note that the Find^ operation on the LR-stack is not needed in the algorithm. 
Therefore, the LR-stack may be easily implemented so as to ensure linear running times for Pop^{bt) 
and Popj^{Bt) with respect to the number dlt and drt of discarded elements, and constant running time 
for PushLB- Then the for loop in steps 4-11 takes running time 0{Yii<t<n{dh + drt + Ft)), where 
Ft is the running time of the Filter procedure for t, that is 0(Si<(<„((i/j + drt) + F). Furthermore, 
each variable bt and Bt is pushed on the L7?-stack at most once (step 8), and thus it is discarded from the 
L7?-stack at most once. Consequently, Si<t<„((iZt + drt) is in 0{n) and the for loop takes 0{n + F) 
total time. Concerning the memory requirements, it is clear that L and R are filled in with elements 
bt,Bt whose cardinality is in 0{n). 

Given the hypothesis that computing the pairs [/?*, Bf] for t e [?^— 1] takes negligible time and space, 
it remains to show that the other computations in steps 1 and 2 take 0{Kn) time and 0{n) additional 
space. To this end, we compute mj^ and Mt, t <E [n — 1], for each permutation Pk in 0{n) time and 
0{n) additional space, as described below. The values computed for each permutation are progressively 
included in the computation of nit, Mt, so as to use a global 0{n) space. 

In |I7]|4J, authors solve a problem called range minimum query (abbreviation: RMQ problem). More 
precisely, they show that, given any array A of n numbers, it is possible to preprocess it in 0{n) time 
so as to answer in 0(1) any query asking for (the position of) the minimum value between two given 
positions qi and q2 in A. This result, closely related to computing the least common ancestor of two 
given nodes in a rooted tree, allows us to compute m^ , Mt for all t in linear time, for each permutation 
Pk- Then we are already done. 

However, we propose here another algorithm, answering a set Q of queries in 0{n+\Q\) time. This 
algorithm is obviously less powerful than the preceding ones, but has at least two advantages. First, 
it is conceptually and algorithmically simpler, allowing the reader to immediately simulate executions. 
Second, it gives another application of LR-stacks, needing this time the implementation of the FindL 
operation. 

Algorithm Compute Inf is given in Algorithm [2] The input is an arbitrary permutation P, signed or 
not, for which we use the same notations as for P^ but without the subscript k. It is quite easy to notice 
that the algorithm works similarly for an array. 

For a pair (qi, 52) with 1 < gi < (72 < "■, denote 

Inf(qi,q2) = m\n{ph \qi<h< 92}. 

Algorithm Computeinf computes Inf{qi, (72) for all pairs in some given set Q, in 0{n + \Q\) time. 
To see this, we first need the following result. 

Claim 2. Let /i G [n\. After the execution of the for loop in step 3 for h, the LR-stack satisfies the 
property i £ Setji{a) ifflnf{i, h) = a. 

Proof of Claim l2j Notice that only allowed operations are performed on the LR-stack, except that 
Pushijfi{ph, h) is not preceded by Popj^{h). This operation would have no effect, since at the beginning 
of the h-t\\ step all the elements in R are akeady less than h. 

Remark that in a i^i?^ -stack, a' is i-blocking for a iff a' > a, and 6' is i?-blocking for b iff 
b' > b. We use induction. For h = Q, at the beginning of the execution the LR-stack is empty, and only 
PushLii{n + 1, 0) is executed, insuring the claim is true. Assuming the claim true for h, let us show it 
is true for h + 1. 



Algorithm 2 The Computeinf algorithm 



Input: Permutation P over [n], Q C {(gi, 52) | 1 < qi < 92 <n} 
Output: Values /«/(gi, (72) for all (gi, (72) G Q 
1: Initialize an empty L~i?~ -Stack 

Po <— " + 1 

for h from to n do 

fo/?^ (p/i) //update Inf{i, h) for alii < h s.t. Inf{i, h — I) > Inf{i, h) 

PushLRiPh, h) lllnf{h, h) = ph 

for all pairs (gi, 92) G Q such that ft, = (72 do 
query{qx,qi) 4- FindL{qi) 



end for 
end for 



Proof of "=>". Let 2 G Setii{a) at the end of the execution of the for loop for h + 1. Two cases may 
appear: 

(i) When a < p/j+i, we deduce that ^(^^^(p/i+i) did not discard a, thus i e Setji{a) at the end of the 
preceding execution of the for loop. By the inductive hypothesis, we deduce that Inf{i, h) = a 
and thus, with a < Ph+i, we have that Inf{i, h + 1) = a. 

(ii) When a = Ph+i, either i has been added by Push]^^{ph^i, h + 1), or some a' > Ph+i existed 
on L before Popj^{ph) such that i e SetR{a') at the end of the execution of the for loop for h. 
In the first case, i — h + I and the conclusion obviously holds. In the second case, the inductive 
hypothesis implies that Inf{i, h) = a' . Consequently, we have that Inf{i, h + 1) = Ph+i = a, 
since a' > Ph+i- 

Proof of "^". By the hypothesis, we assume /n/(i, /i + 1) = a at the end of the execution of the for 
loop for /i + 1. 

(i) When a < Ph+i, we deduce Inf{i, h) ^ a and by the inductive hypothesis we obtain i € Setii{a) 
at the end of the execution of the for loop for h. Since a < Ph+i, the instruction Pop j^{ph+i) 
does not discard a and the conclusion follows. 

(ii) When a = Ph+i, we deduce Inf{i.h) > Ph+i, since Ph+i occurs only once on P. By the 
inductive hypothesis, that means i £ SetR{a') at the end of the execution of the for loop for h, 
with a' — Inf{i, h) > Ph+i- Consequently, Pop]^{ph+i) discards a' and moves i into Setii{ph^i) 
(which is SetR{a)). 

The claim is proved.l 

Claims. Foreach (91,92) G Q, the value que ry{qi,q2) returned by Algorithm Computeinf is equal to 
Inf(qi, (72)- Moreover, the algorithm may be implemented to have 0{n + \Q\) running time. 

Proof of Claim l3| The value query {qi, (72) is computed in step 7 of the for loop, when ft = (72- By 
Claiml2] FindLiqi) = Inf{qi, 92) and we are done. 

Concerning the running time, the difficulty comes from the need to implement the operation Find^. 
However, we benefit here from a very particular case of the Union-Find-Delete context, where no dele- 
tion is performed and the union operations (due to Popj^ and Pushtu) always join sets of consecutive 
elements to obtain another set of consecutive elements. In fT6l, an implementation of the Union-Find 
operations for this particular case is proposed, which realizes each union between two sets in 0(1), and 
m operations Find in 0{n + m). Moreover, the sets (equivalently: the elements of L) as well as their 
elements (equivalently: the elements in R) may be easily chained to simulate the L and R stacks. Then 
each Popi^ operation takes linear time with respect to the number of discarded elements and is in 0{n) 
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Algorithm 3 The Filter algorithm for common intervals 



Input: Pointers R^ (t) , R^ (t) to the first and last element of Setji{t) (possibly equal to nil) 
Output: All common intervals {t..x) of 7^, with fixed t. 
1: if/?^(t) ^ nil then 

x'^ ^ the tai-get of R^ (t); x-'- ^ the tai-get of R^{t) 



X -i— a;^ 

while X < x^ do 

Output the interval {t..x) 

X <— next{x) //or n + 1 ifnext{x) does not exist 

end while 
end if 



time over all h, since the elements on L are elements of P, which are pushed exactly once on L. Fur- 
thermore, PusliLR takes constant time and the |Q| calls of Find l take 0{n+ \Q\) time. The indicated 
running time for our algorithm follows. I 

Proof of Theorem l2] (continued). Algorithm ComputeSup, which computes the maximum values be- 
tween given pairs of positions in P, is clearly similar to Computelnf. Then, m^ and M^ may be obtained 
by appropriately defining the query set Q, and the running time of L7?-Search follows. I 

For each class C among common, nested and conserved intervals, the following sections give the 
settings for b, B and Filter in the L/?-Search algorithm, prove the correction of the resulting variant of 
L/?-Search and discuss subclasses C C C. For the ease of presentation, we assume that R^{a) is a 
pointer to the last element element of SetR{a), for all a G S'L with Setii{a) ^ 0. It is easy to check that 
such a pointer is easily updated during the operations on the L/?-stack. 

Remark 4. According to the preceding definitions, the set Setii{t) is, at the end of step 9 of for^, the 
set denoted Set*j^{t). 

5 Finding common, nested and conserved intervals 
5.1 Common intervals 

In this case, T' is a set of permutations (all elements have a + sign). The M/nMax-profile of V uses in 
this case the basic settings for b and B: 

• b(i) :— nii, for alH, 1 < i < n — 1 

• B{i) := Af„ for all i, 1 < i < n - 1 

• Filter is given in Algorithmic] 

Example 5. Recall that, in Figure[l] V ^ {Pi, P2} with Pi = Idj and P2 = (7 2 1 3 6 4 5). The Filter 
procedure in Algorithm [3] successively outputs the intervals (4. .5), (4. .6) (when t = 4), (3. .6) (when 

t = 3) and (1..2), (1..3), (1..6), (1..7) (when t = 1). 

Theorem 3. Algorithm LR-Search with settings b, B and Filter above solves Common-ISP for K per- 
mutations in 0{Kn + N) time and 0{n) additional space, where N is the number of common intervals 
ofV. 

Proof of Theorem pi By Remark H Algorithm [3] c onsiders each x in Set]^{t). Thus, over all values 
of t. Filter outputs (Ji<t<n-iPairst. By Theorem 111 this is exactly the set of common intervals {t..x) 
of V with bt = t and B^-i = x. Thus all the intervals output by the algorithm are common intervals 
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of V. Conversely, let {t..x) be a common interval of V. Then, for each k, 2 < k < K, and each w, 
t < w < X — 1, we must have t < m^ < w < w + 1 < M^ < x. Consequently, t < m^ < w < 
w + 1 < M^ < X for allw, t < w < X — 1. By definition, bi — b{i) — rrii and Bi = B{i) ~ Mi for 
all i, and thus we have t < b^ and B^ < x for all w. We deduce that 

bt < t < niin{fe^ \t < w < x — 1} 

Bx-i > X > niax{_B^ \t < w < x ~ V\ 

and thus t = bt and x = B^-i- By Theorem [1] we obtain that {t..x) is output by L/?-Search and the 
correction of the algorithm is proved. 

The while loop in steps 4-7 of Filter takes 0(|5ef^(i)|) time, since R^ (t) = top{R). Filter has 
therefore a linear complexity with respect to the size of the output. By Theorem [2] and given that b{i) 
and B{i) are computed in constant time for each i when m^ and Af, are known, the L/?-Search algorithm 
with the abovementioned settings runs in 0{n + N) time. The additional space used by Filter is in 

o(i).l 

5.2 Nested Intervals 

In this case too, T-" is a set of permutations (all elements have + sign) and the M/nMax-profile uses the 
basic settings for b and B: 

• b{i) :— nii, for alH, 1 < ? < n — 1 

• B{i) := M,, for alH, 1 < i < ?i - 1 

• Filter is given in Algorithmffl 

Example 6. Again, in FigurefTlthe LR-Sesn-ch algorithm is applied to T' = {Pi, P2}, where Pi = Id^ 

and P2 ~ (7213645). When t = 4, the Filter procedure in AlgorithmHoutputs intervals (4. .5), (4. .6) 
and sets W[S] ^ 6. When t ^ 3, top{L) = 3 and W[3] ^ 0, so that (3..6) is output and W[2] receives 
the value 6. Now, when t = 2 we have top{L) ^ 2 and the if instruction in step 4 stops the execution 
of Filter. Then W[l\ remains 0. Thus, when i = 1, in step 9 of Filter we have x^ = t + 1 and the 
algorithm outputs (1..2), (1..3) but not (1..6), nor (1..7). 

In this subsection, we say that (x, a; + 1) is a gap if exactly one of x, a; + 1 is on R. 

Theorem 4. Algorithm LR-Search with settings b, B and Filter above solves Nested-ISP for K permu- 
tations in 0{Kn + N) time and 0{n) additional space, where N is the number of nested intervals of 
V. 

By RemarkH Setn.{t) in Algorithmplis the same as Set]^{t). Notice that the value of W[t] computed 
by Filter during its execution f or t + 1 indicates the largest element x in Set*j^^{t + 1) such that {t + 
l..x) is output by Filter (step 13 in Filter for t + 1). However, W[t] may be set to in step 7 of the 
Filter procedure (when executed for t), if we know that its initial value was discarded by Popj^{Bt) in 
the LR-Search algorithm. 

From now on, we focus on the execution of the for loop in LR-Search for t (including Filter). Let 
yt be defined as follows. We use the notations in step 5 of the Filter procedure, and assume thus that 
R^ [t) ^ nil. If x^ =t + l or W[t] ^ 0, then yt is the first element in Set*j^{t) such that the following 
properties hold: 

(1) W[t] < yt, and 

(2) either {yt,yt + 1) is a gap, or yt = x-^. 

Otherwise, i.e. if x^ y^ t + 1 and W[t] ~ 0, yt is set to 0. Let 
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Algorithm 4 The Filter algorithm for nested intervals 



Input: Pointers R^ (t) , R^ (t) to the first and last element of Setji{t) (possibly equal to nil) 
Output: All nested intervals {t..x) of T' with fixed t. 
1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 



i{t — n— 1 then 

Let T/F be a n-size vector filled in with 
end if 

itR'^it) ^ nil then 

a;^ ^ the target of R^ (t); x-^ -(- the tai-get of R^{t) 
if x^ > Ty[f] then 

W[t] -i— //element W[t] has been discarded byPopfi{Bt 

end if 
itx'^ =t + lorW[t]^0 then 



X ^ x^ 



while X < x-^ and (x = t + 1 or x < VK[t] or (x — 1, x) is not a gap) do 
Output the interval {t..x) 

W[t — 1] -i— a; //t passes down tot — 1 its largest x such that {t..x) is output 

X i~ next{x) //or n + 1 ifnext{x) does not exist 

end while 
end if 



17: end if 



Vt -.^ {v e Set'nit) \x^ < V < yt}. 
Now, let us prove that: 

Claim 4. The intervals output by Filter are the intervals (t..x) with x £ Vt- 

Proof of Claim Bl Indeed, when x^ ^ t + 1 and W[t] — 0, the condition in step 9 is false and the 
algorithm returns an empty set of intervals. This is correct, since Vt — %'n\ this case. 

In the contrary case, the condition in step 9 is true. The while loop starts with x ^ x^ and outputs 
all intervals {t..x) with x in Set]i{t) up to a last one {t..xo). If xq = x-^, then xq ~ x-^ = yt, since x-^ 
satisfies properties (1) and (2) in the definition of yt whereas no other element preceding it in Set]i{t) 
does. If xq ^ a;-*-, then xq satisfies the second condition in step 11, whereas the element a:;' — next{x) 
is in Set]^{t) but does not satisfy the second condition in step 11. Then (a;' — l,x') is a gap and thus 
(a;o, Xq + 1) is a gap too. In order to deduce that xq — yt, we only have to show that condition (1) in 
the definition of yt is satisfied by a;o. If, by contradiction, W[t] > xq, then W[t] is on R below xq and 
thus x' < W[t]. This is impossible, since x' does not satisfy the second condition in step 11. 1 

The following claim estabUshes the correction of our algorithm: 

Claim 5. We have x Cz Vt iff(t..x) is a nested interval. 

Proof of Claimpl This proof is by induction on i. Forf = n—l,R^{n—l) ^ niliff PushLii{bt,t + l) 
is called in LR-Seaich with bt = n — 1 and t + 1 = n. Thus Sef^^{n — 1) = {n}. Then, with 
W[n — 1] = 0, we have that V^„_i — {n} and indeed {n — l..n) is a nested interval. 

We now assume the claim is true for i + 1 and show it for t. Then R^ (t) ^ nil in step 4 of Filter 
(otherwise Vt is not defined, and there is no nested interval with left endpoint t) and thus Popj^{bt) in 
step 5 of L/?-Search has been executed with bt = t. Consequently, Vt+i C Setu,{t) at the end of step 5 
in fort of L7?-Search. In step 6, all elements x of R with x < Bt are discarded and, with them, the first 
elements of Vt+i (since the first element of Setfi{t), as well as of Vt+i if it is non-empty, is top{R)). 
Furthermore, the PushLfi{bt , i + 1 ) operation, if executed, pushes t + 1 on the top of R and, in the same 
time, in Set*^{t). 
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When Filter is executed for t, several situations may occur. Notations x^ and x^ concern the 
execution of the for loop in L/?-Search, including Filter, for t. Then x^ = top{R). 

(i) x^ = f + 1 and W[t\ = in step 9 of Filter. Then either there is no x such that [t + l..a;) is 
output, or such an x exists but is removed from R by Popj^{Bt). According to the definition, yt 
is the first element of Set]^{t) such that either (j/i, j/t + 1) is a gap, or yt = x-^. In both cases, 
since x^ = t + 1 we have that i + 1 is in Vt- Also, {t..t + 1) is nested since it is a common 
interval (by TheoremfT) and has cardinality 2. By the definition of nested intervals, all intervals 
{t..x) with successive values of x, a; > t + 2, will be nested, thus x E Vt implies that {t..x) 
is nested. Conversely, assume by contradiction that some x > t + 2 exists such that {t..x) is 
nested but x ^ Vt, and take x as small as possible with these properties. Since {t..x) is nested, 
{t..x) is a common interval thus, by Theorems fl] and [3] and given the settings of b and B, we have 
X G Set^f^it), thus X < x-^. As {t..x) is nested, either (t + l..x) or {t..x — 1) is nested. The 
former does not hold (because of W[t] = 0), thus {t..x — 1) is nested. Since x was the smallest 
with the indicated properties, x — 1 G Vt- But then, by the definition of yt, we have x E Vt too, a 
contradiction. 

(ii) x^ ^t + landW[t]>0 in step 9 of Filter. Then 

Vt = {t + l}U Vt+i -{veVt+i\v<t+l}UU 

where U is the set of all consecutive elements W[t] + 1, W[t] + 2, . . . ,z C Set]^{t) such that 
either z = x-^ or {z, z + 1) is a gap. Equivalently, Vt gets all the elements in Vt+i except those 
smaller than t + 1 (because of Pop j^), as well as t + 1 (pushed by PusIilh) and all the consecutive 
elements that are possibly added at the end of Vt+i during the Popj^ipt) operation, with bt = t. 
It is easy to see that the last element of U satisfies the conditions (1) and (2) in the definition of 
yt, and no one before it in Setj^{t) does. Thus z ^ yt- 

By contradiction, assume some x E Vt exists such that {t..x) is not nested. Assume x is the 
smallest with these properties. Then x > W[t], otherwise x = t + 1 or a: e Vt+i and then 
(t + l..x) is nested by definition and the inductive hypothesis, insuring that {t..x) is nested. Now, 
X — 1 E Vt (since there is no gap in Vt beyond W[t]) thus, by the minimality of x, {t..x — 1) is 
nested. But then {t..x) is nested, a contradiction. Conversely, let {t..x) be a nested interval, and let 
us show that x G Vt. Once again, assume this is not true and let x be the smallest counter-example. 
Since {t..x) is nested, {t..x) is a common interval and, by Theorems n] and pi x G Set\^{t) thus 
X > t + 1. Now, we must have x > W[t], otherwise x = t + loixE Vt+i thus x E Vt, ^ 
contradiction. Finally, since {t..x) is nested we have two cases. Either [t + 1, x) is nested, and 
then by the inductive hypothesis x E Vt+i thus x E Vt,a contradiction. Or (t, a; — 1) is nested and 
thus x — 1 € Vt by the minimality of t, thus x e Vt by the definition of Vt, another contradiction. 

(iii) x^ 7^ t + 1 and W[t] > in step 9 of Filter. Then 

Vt - Vt+i -{vEVt+i\v<t+l}UU 
with U defined as previously done, and the proof follows similarly. 

(iv) x^ ^ t + 1 and W[t] — in step 9 of Filter. Then Vt = and we must show there is no nested 
interval {t..x). If, by contradiction, such an x exists, then assume x is taken to be the smallest 
one. Then {t..x — 1) is not nested, by the minimality of x. Thus {t + 1, x) is nested, and thus 
X E Vf+i. But then x E Vt unless x is removed by Popj^{Bt), with Bt = t + 1. However, this is 
impossible, since a; > i + 1.1 

Remark 5. According to the preceding claim, W[t] = yt+i- Moreover, assume that every element r 
in R has an associated pointer R3{r) on the first element w E R larger than or equal to r and such that 
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either {w, u; + 1) is a gap or w is the bottom of R. Then yt may be computed in constant time in each 
of the cases (i)-(iv) of the proof, using W[t] (i.e. yt+i) and _R^(ro), where tq is the target of R^ (t) at 
the end of fort+i, if R^ (t) =^ nil. 

Proof of Theorem SI The algorithm correction is proved by Claims [4] and Is] The 0{Kn + N) running 
time of the algorithm is due to Theorem [2] and to the linearity of Filter with respect to \Vt\, which is 
clear f or t < n — 1 since Filter stops when the first element not in Vt is found. The 0{n) complexity 
when i = n — 1 does not change the overall running time of the algorithm. I 



5.3 Conserved intervals 

Here, each Pj, is a signed permutation with p\ = 1 and pj^ = n, both positive. The definitions of b and 
B define a M/nMax -profile of V adapted to the specific needs of conserved intervals. For each i with 

1 < i < n — 1, let u^ :~ m^ if m^ — i, and u^ := rn^ — 1 otherwise. Similarly, let v'^ := Mj^ if 
M^ — i + 1^ and wf := Mj' + 1 otherwise. Let: 

• b{i) := inm{u^ | 2 < fc < K}, for all i, 1 < i < n - 1 

• B{i) := max{i;f | 2 < fc < iiT}, for alH, 1 < i < n - 1 

• Filter is given in Algorithm pi where: R^ {t)* is a pointer to the first element x in Setji{t) such 
that t and x have the same sign in all permutations in V (such elements x are chained together 
inside R); and Position{t. t + 1) returns true iff, for each fc, either t is positive in Pk and Pf^^{t) < 
Pk^{t + 1), or t is negative in Pk and P^T^t) > P^^it + 1). 

Remark 6. Note that {t..x) isa conserved interval of V iff it has the following properties: 

(1) it is a common interval of V 

(2) it is delimited by t and a; on P^:, for all k E [K] 

(3) t and x are both positive or both negative in each Pj., for all k E [K] 

(4) for each k E [K], either t is positive in P^ and Pj7^{t) < P/^^ix), or t is negative in Pj. and 

Conditions (1) and (2) are easily handled by defining the bounding functions b and B as indicated. 
However, conditions (3) and (4) need a preprocessing of the permutations in V, in order to: (Task 1) 
identify and chain together inside R the elements x in the same equivalence class with respect to the 
relation "x and x' have the same sign in all permutations", allowing us to deal with (3); and (Task 2) 
compute the boolean function PositionQ defined above, which will insure that (4) holds. These tasks 
are done in 0{Kn) time and 0{n) additional space as follows: 

Task 1. Consider the matrix M whose row vectors are the vectors sigup for fc E {2, . . . , K}, 
and perform a radix sort on the columns of this matrix (which correspond to the elements t of the 
permutations). Group together all elements t that have the same column vector, i.e. the same sign in 
all permutations. For each group s, the elements of the group that are pushed on R are progressively 
chained together immediately after PushLfi{bt, t + 1) (step 8 in LR-Seaich), and the pointer First{s) to 
the first element of the chain is updated as needed. Then, for each t, R^ (t)* is defined as First{s), where 
s is the group of i. All these operations are done in 0{Kn) time and 0{n) additional space, assuming 
that radix sort does not really create the indicated matrix, but rather manipulates column numbers and 
checks the values directly on the vectors signp (which are not modified). 

Task 2. The relative positions of t and t + 1 in each Pj. are checked in 0(1) using the functions 



Pf. ^(), and thus Poj'/f/on(t,i+l) is computed in 0{K) time and 0(1) additional space, for each f E 



With this supplementary information, the Filter procedure is presented in Algorithm|5] 
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Algorithm 5 The Filter algorithm for conserved intervals 



Input: Pointers R^ (t)* to the first element in Setji{t) having the same sign as t in all permutations Pk, 

and R'^{t) to the last element in Setji{t) (they are possibly nil) 
Output: All conserved intervals {t..x) of V with fixed t. 

1: iiR^{t)* T^mHhen 

v^ -f- the tai-get oiR^ {t)*; x^ <- the target of R^{t) 



a: <— x^ 



if Position{t, t + 1) tlien 
wliile a; < a;^ do 
6: Output the interval {t..x) 

7: X -^ the element immediately following x in its chain //or n + 1 if it does not exist 

8: end wliile 

9: end if 
10: end if 



Example 7. In Figure l2J conserved intervals are computed for V — {Pi,P2}, where Pi = Idj and 
P2 — (13 2 6 4 5 7) with positive 1, 6 and 7, and negative 2, 3, 4, 5. Note, for, instance, that rng = 4 
but bg — 3, by the definition of bt, indicating that 4 cannot be the delimiter of a conserved interval 
containing 6 and 7. Similar remarks are vaUd for 65, i?3, 63 and Bi. During the execution of L/J-Search, 
Filter does not output (4. .5) since Position{A,5) — false, but outputs (2. .3) and (1..7) for which 
variable Position{) is true and the signs are compatible. These are the conserved intervals of T'. 

Theorem 5. Algorithm LR-Search with settings b, B and Filter above solves Conserved-ISP for K 
signed permutations in 0{Kn + N) time and 0{n) additional space, where N is the number of con- 
served intervals ofV. 

We start by showing that: 

Claim 6. An interval {t..x) with the properties (1), (2) in Remarkpjalso satisfies property (4) iff 
Position{t, i + 1) is true. 

Proof of Claim [6j The "=>" part is obviously true. For the "<^" part, assume by contradiction and 
without any loss of generality that k and x, with x > t, exist such that t is positive and Pi7^{x) < 
P^^{t). Then, the hypothesis that Position[t,t + 1) is true insures that P^^{t) < P^^^{t + 1), thus 
t + 1 does not belong to the interval delimited by t and x. Consequently, condition (1) in Remark l6] is 
contradicted. The reasoning is similar if t is negative. I 

Proof of Theorem I5J We show that {t..x) is output by L7?-Search with the given parameters iff it is a 
conserved interval. Recall that, by RemarkH Setj^{t) in AlgorithmBlis the same as Set*'j^{t). We assume 
without loss of generality that t is positive. 

Proof of "=>:" According to Theoremfll the set of intervals computed by the algorithm L/?-Search 
is the set of common intervals (ignoring the signs) {t..x) of V with bt — t and B^-i — x. Then, for 
each output {t..x) and each k, the corresponding interval on Pk has property (1) in Remarkl6] To check 
property (2), notice that by TheoremfTlwe have t — bt and x — B^-i, thus according to the conditions 
on^()andB(): 

mt < t = bt = b{t) < nit and M^^i > x ~ B^-i = B{x — 1) > M^-i 

We deduce that bt — mt — t and B^-i = M^-i — x. By the second part of Theoremfl] we know that 
bt < biu and B^^i > Byj for ?i\\w,t<w<x—\. Assume by contradiction that t is not a delimiter 
of the interval of Pj. made of t,t + 1, ... ,x. Then, there is some w,t + \<w<x — 1 such that t is 
between w and w + 1 on Pk- Then m.u, < m^ = t < w and thus fe^ = b{w) < ra^ — 1 < i — 1. But 
then bt = t>t— 1> bw, a contradiction. The reasoning is similar for x. Property (2) is proved. 
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+ [1,4] [2,3] [1,7] [4,5] [3,6] . [3,7] . 

LR LR LR LR LR LR 



LR 



2 — >3 1 — >7 4 >5 3 — >6 3 

1^7 3^6 7 

7 



PoPi(l) Popt(2) Popi(l) Popt(4) Popt(3) Popi(3) 

P0Pr(4) P0Pr{3) P0Pr(7) PoPr(5) PoPr(6) PoPr(7) 

Push 1^1^(2,3) Pus/i^R (4,5) Pus/i ^^(3,6) Push i^Ci.l) 



Figure 2: The L/?-Search algorithm for conserved intervals when P2 = (13 2 6 4 5 7) with the signs 
indicated on the figure. The chains in R may be easily deduced using the signs. Bounds bt and Bt are 
in bold font whenever they are not equal to nit and respectively Mt- 

Property (3) is insured by the interpretation of R^ (t)* and R'^{t), as well as by steps 1-3 and 7 in 
Filter. Claim [6] and step 4 in Filter guarantee that the property (4) holds. Thus {t..x) is a conserved 
interval. 

Proof of "4=:" We have bt = b{t) = nit ^ t and B^-i = B{x - 1) = Af^-i = x, by the definition 
of b, B and since a conserved interval is a common interval. By Theoremll] we deduce that x e Set*j^{t). 
Furthermore, we use properties (l)-(4) in Remark |6] and show that no interval with these properties is 
forgot by Filter. By contradiction, if this was the case, then x would be eliminated by the condition in 
step 4 of Filter. But then, by Claiml6]the interval {t..x) cannot satisfy property (4) in Remarkl6] This 
contradicts the assumption that {t..x) is conserved. 

It is easy to see that bt and Bt may be computed in 0{Kn) time and 0{n) additional space using the 
values mt,Mt (and avoiding to store all these values), thus we may apply Theoremll Filter clearly has 
running time proportional to the number of output intervals, since the while loop in step 5 is executed 
only when the condition in step 4 is true. Moreover, the while loop has running time proportional to the 
number of output intervals. Theorem0finishes the proof.l 



6 Finding subclasses of common, nested and conserved intervals 
6.1 Irreducible common intervals, and other common intervals 

Irreducible common intervals have been defined in | [T2| as follows. Let G be the graph whose vertices 
are all the common intervals of V, and whose edges are the pairs of non-disjoint common intervals. 
Then a common interval / is reducible if it is the set union of some of its proper sub-intervals that are 
common intervals and induce together a connected subgraph of G. Otherwise, / is irreducible. 

Consider now the total order on the set of common intervals of P given by {ti..xi) < {t2.-X2) iff 
either ti > t2, or ti — ^2 and xi < X2- For each w with 1 < ?« < n — 1, let Small (w) denote the 
smallest, with respect to this order, common interval of V containing w and w + 1. It is shown in lil2J 
that: 

Claim 7. The set of irreducible intervals ofV is the set {Small (w) \ 1 < w < n — 1}. 

Then lfT2l proposes an algorithm to solve IrreducibleCommon-ISP in linear time. Another linear 
time algorithm is obtained by appropriately filtering the results of our LR-Search algorithm, as shown 
below. 
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Claim 8. Let w be an integer such that \ < w < n — 1. The interval Small (w) is the common interval 
{t..x) ofV such that the couple (n — t, x) of integers is minimum with respect to the lexicographic order 
with the property t < w < x. 

Proof of Claiin|8j The condition involving the lexicographic order means that we first maximize t, and 
then minimize x. By contradiction, assume that the interval {t..x) is not Small (w). Now, Small (w) 
is minimum according to the total order. We must then have either t < min{Small {w)) (the smallest 
element in Small (w)), or i = min{Small (w)) and x > maxiSmall {w)). But this contradicts the choice 
of t and X.I 

The L/?-Search algorithm we propose here uses the same MmMax-profile for V as common intervals, 
that is b{i) := rrii and B{i) :— Mi. However, the Filter algorithm is in this case a refinement of the 
initial Filter procedure in Algorithm [3] Algorithm l6] uses a stack S to store the values w for which 
Small [w) has not been found yet, in increasing order from top to bottom of S. For a given w, the first 
interval {t..x) such that t < w < a; found by LTJ-Search is Small (w). However, such a value x may be 
located in Set*'p,{t) after a sequence of useless values x' . At each moment of the execution of Filter, say 
that a value x on i? is untrusty if at least one interval {t' ..x) has been already output (meaning that x 
has possibly become useless), and trusty otherwise. Call a strip any maximal sequence of consecutive 
integers on R that are all untrusty. 

In order to insure the best running time for our algorithm, we consider that each element a; on i? but 
the bottom of R has, in addition to its successor next{x) on R, another successor denoted nextt{x): 

• if a; is trusty, then its successor nextt{x) is next{x); 

• if X is untrusty and is the head (i.e. the first) element in its strip, then nextt{x) is the first element 
following the strip (if such an element exists) 

• if X is untrusty and it is not the head element in its strip, then nextt{x) — nextt(x'), where x' is 
the head element in the strip of x. 

Remark 7. Note that the state (trusty or untrusty) of an element may be easily computed in Filter. 
Moreover, the successor nextt{x) of an element x is computed in 0(1) when x is pushed on R by 
Push^R- Furthermore, when an element x becomes untrusty (step 10 in Filter), the successors nexttQ 
change in the strip before and possibly in the strip after x on R. We call Update(nexff()) the procedure 
performing this update in Algorithm l6] (step 10). We do not give its details here, but discuss it in the 
proof of Theoreml6] 

Example 8. On the example in Figure [Tl the L/?-Search algorithm with the Filter procedure in Algo- 
rithml6]first outputs (when i = 4) (4. .5) and (4. .6) (which are Small (4) and Small (5)). The values 5 
and 6 become untrusty (step 10). When t — 3, the interval (3. .6) is output, which is Small (3). Finally, 
with t = 1, intervals (1..2) and (1..3) are successively output (they are Small (1) and Small (2)). The 
stack S still contains 6 and the next value in SetR{l) is x = 6. Thus top{S) = x = 6 in step 8, and 
X is untrusty. In step 16, x is updated to the value nextt{Q) which is 7. The interval (1..7), which is 
Small (6), is output. 

Theorem 6. Algorithm LR-Search with settings b,B as for common intervals and the more restrictive 
Filter procedure in Algorithmplsolves IrreducibleCommon-ISP for K permutations in 0{Kn) time and 
0{n) additional space. 

Proof of Theorem |6J The decreasing order of t in the for loop of the L7?-Search algorithm and the 
increasing values of the elements in Set]^{t) show that the intervals (t..x) are considered according to 
the lexicographic order required by ClaimlS] Values w are pushed on S as soon as w is considered (step 
4 in Filter) and are discarded from S iff the interval with maximum t and minimum x containing w and 
w + 1 is output (steps 10 and 12). By Claimls] this interval is Small (w). Thus, all the intervals output 
by the algorithm are irreducible. 
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Algorithm 6 The Filter algorithm for irreducible common intervals 



itt = n— 1 then 

Let S be an empty stack 
end if 

Push t on S //Small (t) has not been found yet 

itR^{t) ^ nil then 

x'^ ^ the tai-get of R^ (t); x-^ -s- the tai-get of R^{t) 
X •<— x^ 

while X < x-^ and S is not empty and {top{S) < a; or a; is untrusty) do 
if fop (5) < xthen 

Output the interval {t..x); Update(nexff()) //x becomes untrusty 

while S is not empty and top{S) < a; do 

Pop top{S) from S // fortop{S), interval Small {top{S)) has been output 

end while 

X -!— next{x) //or n + 1 ifnext{x) does not exist 

else 

X -!— nextt{x) //orn + 1 ifnexttix) does not exist 

end if 
end while 
end if 



Conversely, assume by contradiction that some interval Small (w) is not output by the algorithm, 
and let w be the smallest such value. Let Small (w) = {t..x*) and consider the execution of Filter for 
t. Since {t..x*) is not output by the algorithm, the execution of the while loop in step 8 either (i) misses 
X* by skipping it in step 16, or (ii) stops before x* is reached. Notice that x* € Set^f^{t) (by Theoremllb. 
Let x' be the largest element smallest than x* for which the condition in step 8 of Filter is tested, and 
consider the state of the stack S immediately after this test. 

Let us show that top{S) = w. By contradiction, if we assume u :— top{S) < w then the mini- 
mality of w insures that Small (u) is output by the algorithm LR-Seaich with the given settings. Then 
Small (u) = (to..a;o) with Iq < t, since u is still on S. Moreover, u > t since Pushj^ji has pushed on 
R only elements t' + 1 with t' > t. Now, we cannot have to < t since then to<t<u<'w<x* 
and thus {t..x*) is smaller than Small (u) and contains both u and u + 1, a contradiction. We thus have 
to = t. Now, {t..xo) is certainly output by Filter, and this has not been done yet when x' is considered 
(otherwise, u would have been discarded from S). We deduce that x' does not stop the execution of 
the while loop, thus m < a;' or x' is untrusty. The former one would contradict the maximality of x', 
because of step 14 which brings into step 8 a value larger than x' . The latter one implies that the next 
trusty value is larger than x* (since x* is skipped because of x' <— nextt{x')), and thus xq > x* and 
{t..x*) < Small (u) contradicts the minimality of Small (u). Thus top{S) = w. 

Consider now the cases (i) and (ii) before. 

(i) In this case, x' is the value for which step 16 has been executed. Then x' and x* are untrusty (by 
the definition of nexttQ), x* > x' (otherwise x* is not skipped) and the values in Set]^{t) between 
x' and X* are consecutive and untrusty (they belong to the same strip). Moreover, step 16 has 
been executed, so that the condition in step 9 of Filter is not verified. Thus top{S) > x'. Since 
top{S) — w and w < a;* (recall that {t..x*) contains w and uj + 1) we deduce that x* > w > x'. 
Then w is between x' and a:* — 1 on the strip, and ui + 1 is between x' + 1 and x* on the strip. 
Consecutively, i(j + 1 is untrusty, thus an interval {t'..w + l) withf' > f has been already output by 
Filter. But this interval is smaller that {t..x*) and should therefore be Small {w), a contradiction. 
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(ii) In this case, x' satisfies x' < x* and x' < top{S) and x' is trusty. Now, x' — 1 is not on S, 
since w ~ top{S) and elements in S are in increasing order from top to bottom. Thus, by the 
minimahty of w. Small {x' — 1) is output by the algorithm. Let Small [x' — 1) = [t" ..x") and 
notice that x" > x' (x' must be contained in Small {x' — 1) but x' has not become untrusty) and 
t" > t (otherwise, x' — 1 should be on S). Then {t..x') n [t" ..x") = [t" ..x') and it is a common 
interval containing both x' ~ \ and x' and which is smaller than Small (x' — 1), a contradiction. 

The correction of the algorithm is proved. We discuss now the running time of the algorithm. 
Leaving apart temporarily the update of nextt{), the running time of Filter is proportional to vt, where 
Vt is the number of elements w for which Small [w) has smallest value t (these elements are discarded 
in step 12). Over all the executions of Filter, we obtain 0{v), where v — I]i<t<„wt. Now, v is in 0{n), 
since it counts the total number of elements w. 

Weneednow to show that the update of nexff() instep lOof Filter may be done in 0(1). For this, we 
see each strip s as a set T{s) which has a representative element rT{s) such that nextt{rT(s)) gives the 
successor of all the elements in the strip s (as insured by the definition of nexttQ). Then it is sufficient 
to update nextt{rT(S)) in order to update nextt{) for all the elements in s. Now, notice that strips may be 
changed in two ways: 

(1) the instruction Popj^ in L/?-Search may perform deletions from some strip s. However, we keep 
the deleted elements in the set T{s) and thus the representative ?'t(s)> ^s well as the set T{s), are 
unchanged by these deletions (although the strip itself is reduced). 

(2) in step 10 of Filter, the element x becomes untrusty and then the strip immediately preceding x (if 
any), the strip formed by x alone, and the strip immediately following x (if any) may concatenate 
(altogether or only two of them, which have consecutive elements). Concatenations of strips 
imply unions of the corresponding sets. 

In conclusion, updating the strings and being able to find the representative of each of them places us 
in the context of a Union-Find structure. The sets are the sets T{s), containing the elements of s as well 
as all the elements previously in s and discarded from R by Popj^. The unions between sets are given by 
the concatenations between strips, whereas the find operation for an element y seeks the representative 
element of the strip containing y. Again, we are in the particular case where unions always involve sets 
of consecutive integers (this is the case when two strips are concatenated). Thus, according to the result 
in lfT6l . the implementation of unions and finds may be done in time linear with respect to the number 
of union and find operations. 

Consequently, the operation Update(nexff()) in step 10 needs to concatenate if necessary two or 
three of the strips indicated in (2) and to update nextt{) for their representative elements. Over all the 
executions of Filter, these operations are done in 0{n) and the running time of the algorithm is proved. I 

Other classes of common intervals 

In [ 13 1, an algorithm is proposed for finding, in K signed permutations, all common intervals whose 
elements have the same sign inside each permutation. To solve this case, our LR-Sesuch algorithm needs 
to (1) preprocess V to compute, for each t and k, the minimum x^, > t, such that t and x^. have different 
signs in P^; and (2) stop to output intervals in the Filter procedure (step 4 in Algorithm [3]l as soon as 
the first X with x > min{a;^. \2 < k < K} is found. 

For a fixed k, task (1) is easily done in 0{n) by considering the elements of Pk in decreasing order, 
and remembering the sign changes. The global time required for this task is thus in 0{Kn), whereas 
the additional space may be limited to 0{n) by computing the min() values above progressively. 

6.2 Maximal nested intervals 

In IIT4II . authors define a nested interval / to be maximal if it is not included in a nested interval of 
size |/| + 1. In fSl, an efficient algorithm is proposed to solve MaximalNested-ISP for K = 2. Not 
surprisingly, LR-Seaich works in this case too, with an appropriate filtering algorithm. 
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To this end, note that: 
Claim 9. A nested interval {t..x) is maximal iff it satisfies the two following conditions: 

(a) (x, X + 1) is a gap with x <E Vt, or x — yt 
(6) bt-i <t~ 1 or Bt-i > X. 

Proof of Claim [9} B y definition, {t..x) is maximal iff it is nested but neither {t..x + l) nor (t~ l..x) are 
nested. By Claimp] this is equivalent to x £ Vt but x + 1 ^ Vt and x ^ Vt-i (or Vt-i is not defined). 
Now, x G Vt and x + 1 ^ Vt simultaneously hold iff property (a) in the claim is satisfied. 

Moreover, x G Vt and x ^ Vt-i (or Vt-i is not defined) iff, at the end of step 9 of fort_i, either 
t — 1 is not on top of L or x has been removed from R. Equivalently, i — 1 > ht-i or x < Bt_i. The 
claim is proved.l 

Example 9. It is easy to see that in FigurefTl only intervals (3. .6) and (1..3) satisfy these conditions. 
They are indeed the only maximal nested intervals of V. 

Then we have: 

Theorem 7. Algorithm LR-Search with settings b,B as for nested intervals and an appropriate Filter 
algorithm solves MaximalNested-ISP for K permutations in 0{Kn + N) time and 0{n) additional 
space, where N is the number of maximal nested intervals ofV. 

Proof of Theorem l7j According to Claim l9] Filter should output only the intervals that satisfy condi- 
tions (a) and (5). The latter condition is easy to test. The former one needs to find each gap, as well as 
yt, in 0{\). For this, it is sufficient to compute and store, for each r E R, the pointer R^{r) defined in 
RemarkB] Values R3{t + 1) must be initialized immediately after PushLii{bt,t + 1), in step for^ of 
LR-Seaich. They do not need to be updated. 

Then it is sufficient to modify Filter in Algorithm |4] by replacing steps 10-15 with the following 
ones: 

X ■<— R^{x^); Compute yt 

while X < yt and (6t_i < f — 1 or i?t_i > x) do 
Output the interval {t..x) 

W[t — 1] <— a; //t passes down tot — 1 its largest x such that {t..x) is output 

X <— R^ {next{x)) //orx <— n + 1 ifnext(x) does not exist 



end while 



Computing yt in 0(1) is possible according to Remark [s] The other modifications aim at precisely 
selecting the elements with properties (a) and (b) in ClaimloF Notice that when a first value x not satis- 
fying properties (a) and (b) is found in step 11, it is clear that no other subsequent value x' (necessarily 
x' > x) will satisfy properties (a) and (6). The resulting Filter procedure then finds all the maximal 
nested intervals, by Claim l9j runs in global time proportional to the number of output intervals, and 
globally uses 0{n) additional space.l 

6.3 Irreducible conserved intervals 

In Q, authors define a conserved interval to be irreducible if it is not the union of smaller conserved 
intervals. They also give an efficient algorithm to solve IrreducibleConserved-ISP. Such an algorithm 
may also be obtained using L/?-Search and the following easy result: 

Claim 10. Let (t..x) be a conserved interval ofV. Then {t..x) is irreducible iff 



min{/i \t < h and {t..h) is a conserved interval of T'}. 
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Proof of ClaimllOJ In |6 1 it is shown that two different irreducible intervals are either disjoint, or nested 
with different endpoints or else overlapping on one element. The conclusion follows. I 

We deduce: 

Theorem 8. Algorithm LR-Search, with settings b, B as for conserved intervals and a simplified Filter 
procedure, solves IrreducibleConserved-ISP for K permutations in 0{Kn) time and 0{n) additional 
space. 

Proof of TheoremlH} By Clai m[TO] it is sufficient to replace the while loop in the Filter procedure with 
an instruction that outputs {t..x ' ).l 

7 Conclusion 

The L/?-stack we introduced in this paper is a simple data structure, of which we noted at least two 
advantages: it is powerful (we had two applications of it in this paper), and it is algorithmically efficient, 
since it makes use of the efficiency reached by the Union-Find-Delete algorithms. 

Using L/?-stacks, our algorithmic framework L/?-Search succeeds in proposing a unique approach 
for dealing with common intervals and their subclasses. The computation of the interval candidates is 
driven by the M/nMax-profile and the bounding functions, whose role is to guarantee that all interval 
candidates satisfy the content-related constraints. Afterwards, the Filter procedure chooses between the 
candidates those that satisfy the supplementary constraints defining a precise subclass. 

All the algorithms resulting from this approach are as efficient as possible. They allowed us to prove 
the power and the flexibility of our approach. Among them, the algorithms searching for nested and 
maximal nested intervals of K permutations, with K > 2, solve previously unsolved problems. 
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